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ABSTRACT 

Using the Submillimeter Array (SMA), we have imaged the 1.3 milhmeter continuum emission 
at the center of the massive star- forming regions NGC 6334 I and I(N). In both regions, the SMA 
observations resolve the emission into multiple millimeter sources, with most of the sources clustered 
into areas only 10,000 AU in diameter. Toward NGC 6334 I, we find four compact sources: the 
two brightest (I-SMAl and I-SMA2) are associated with previously-known ammonia cores; I-SMA3 is 
coincident with the peak of the compact Hii region (NGC 6334 F); and I-SMA4 is a newly-discovered 
object. While I-SMA3 exhibits a mixture of free-free and dust emission, the rest of the objects are 
dust cores. Toward NGC 6334 I(N), seven compact dust cores are found, one of which is associated 
with a faint centimeter source. With the exception of I-SMA3, none of the millimeter sources have 
infrared counterparts in Spitzer Space Telescope 3-8 iim images. Using a simple physical model for the 
dust continuum emission, the estimated mass of interstellar material toward each of these compact 
objects is in the range of 3 to 66 Mq. The total mass in the compact objects appears to be similar 
in I and I(N). The small size of these groups of sources suggest that these objects are proto- Trapezia 
forming in the centers of clusters of low to intermediate-mass stars. 

Subject headings: stars: formation — infrared: stars — ISM: individual (NGC 6334 I) — ISM: indi- 
vidual (NGC 6334 I(N)) — techniques: interferometric — submillimeter 



1. INTRODUCTION 

The formation process of massive stars continues to be 
a poorly-understood phenomenon in astrophysics. The 
most fundamental clues to the origin of OB stars are 
their multiplicity and their common association with 
high column densities of molecular gas and dust. A 
unique feature of OB stars is that they are often found in 
non-hierarchical (and consequently) non-stable systems 
in the center of clusters (Sharplcss 19 5-^ . A nearby 
example of this phenomenon is the Orion Trapezium 
with projected stellar separations of 4000 to 10000 AU, 
which in turn is found in the center of a rich cluster 
of low mass stars. A long-standing question is whether 
these multiple systems are the direct result of the for- 
mation process, thus providing a clue to the process 
of high mass star formation, or the result of dynami- 
cal evolution in the centers of young clusters where the 
most massive stars move toward the center by eject- 
ing lower mass stars outward(Bonnell & Davies 199^. 
While recent (sub)millimeter studies have identified good 
candidates for massive protoclusters on scales of sev- 
eral parsecs, the identification of proto- Trapezia is a 
more difficult prospect. Due to the complex nature of 
high mass star-forming regions, the high extinction typi- 
cally observed toward massive protostars, and their typ- 
ically large distance (> 1 kpc), high-resolution imag- 
ing at (sub)millimeter wavelengths is required to resolve 
one protostar from another. Recent space-based near- 
infrared (NIR) imaging has succeeded in identifying a 
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5600 AU diameter cluster of five protq -OB stars com- 
prising W3 IRS 5 l)Megeath et al J 12005(1 . However, be- 
cause massive protostars form in the deeply-embedded 
cores of molecular clouds, dust extinction may obscure 
a significant fraction of them in the infrared. In these 
cases, millimeter continuum emission from dust pro- 
vides one of the few alternative tracers of protostars be- 
cause it remains optically thin at high column densities 
(A^H < 1025cm-2). 

The recent commissioning of the Submillimeter Array 
(SMA) on Mauna Kea, Hawaii has expanded the range of 
millimeter interferometry to higher frequencies and lower 
declinations. Now within reach is NGC 633 4 (6 = —35°) , 
a luminous and relatively nearby (1.7 kpc: lNeckelll978f) 
molecular cloud/Hii region complex containing sev- 
eral concentrations of massive star formation at various 
stages of evolution IjStraw fc Hvlandlll989() . The north- 
eastern end appears to be the youngest and co ntains 
the ra dio source "F" (G351.42+0.64; Rodrfgucz lit alJ 
p^8^) which is associated with IRAS 17175-354 4. The 
earliest far-infra red images of this region ( Emerson et alJ 
[T973t iMcBreen et all Il979l) identified the emission as 
source "I" . An additional componen t "I(N)" was first de- 
tected at Imm l)Cheung et al.lll978l) and later at 400/xm 
(|Gezarilll982t) . Further observations have demonstrated 
that although these two cores have comparable mass, 
NGC 6334 I dominates the combined bolomctric lumi- 
nosity of 2.6 X 10^ Lq l|Sa^dell 2000) . In the NIR, an 
embedded cluster of st ars has been detect ed in the cen- 
tral 2' of NGC 6334 I l|Tapia et all 1199^ . In contrast, 
the only NIR emission detected toward I(N) are H2 knots 
that are mo st likely associated with outflo w activity (see 
Figure 3 of IMegeath fc Tieftrund l|1999(l ). In this pa- 
per, we present the first 1.3 millimeter interferometry of 
NGC 6334 I and I(N). Despite their strikingly different 
appearance in the NIR, we find that both regions contain 
a similar cluster of compact dust continuum cores. A de- 
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tailed treatment of the millimeter line emission accompa- 
nying these objects will be the subject of a forthcoming 
paper (Hunter et al., in preparation). 

2. OBSERVATIONS 

2.1. Suhmillimeter Array (SMA) 

The SMA^ observations were made with six antennas 
in both the compact configuration (2004 May), and ex- 
tended configuration (2005 May). Two pointings were 
observed: NGC 6334 I at 17''20'°53".44, -35°47'02".2 
and NGC 6334 I(N) at 17'^20'"54".63, -35°45'08".5 
(J2000). Unprojected baseline lengths ranged from 22m 
to 226m. The SMA receiver s are double sid eband mix- 
ers with 2 GHz bandwidth l|Blundel]l l200l . The cen- 
ter frequencies were 217.6 in LSB and 227.6 GHz in 
USB. The correlator was configured for uniform resolu- 
tion (0.81 MHz per channel). Typical system tempera- 
tures were 200K. The gain calibrators were NRAO 530 
(23° distant) and J1924-292 (27° distant) and the band- 
pass calibrators were Uranus and 3C279. The data were 
calibrated in Miriad, then exported to AIPS where the 
line and continuum emission were separated with the 
task UVLSF. Self-calibration was performed on the con- 
tinuum data, and solutions were transferred to the line 
data. Flux calibration is based on SMA flux monitor- 
ing of the observed quasar and is estimated to be accu- 
rate to within 20%. The estimated accuracy of the abso- 
lute coordinates is 0.4". After combining the calibrated 
LSB and USB continuum uv-data, the \a rms noise level 
achieved in the continuum images is 7 mJy beam~^. The 
resulting synthesized beam is 2'.'1 x 1'.'2 and the primary 
beam is ~ 56". 

2.2. Very Large Array (VLA) 

Archival 3.6 cm data from the NRAO^ Very Large Ar- 
ray (VLA) were calibrated and imaged in AIPS. The ob- 
servation date was May 7, 1990 (A-configuration). The 
fl-ux calibrator was 3C286, and the gain calibrator was 
NRAO530. The synthesized beam of the 3.6 cm contin- 
uum image is 0.90" x 0.43" at (P.A.=-18°) and the rms 
sensitivity is 1 mJy bcam^^. 

2.3. Spitzer Space Telescope 

Both SMA field s were observed with the IRAC camera 
ijFazio et alJl200^ onboard the Spitzer Space Telescope. 
The total integration time was 666 seconds. The high 
dynamic range mode was used to obtain 0.4 and 10.4 
second integrations. The data were reduced using the 
BCD data from the Spitzer Science Center version 11.4 
pipehne. When constructing images, the 0.4 second data 
were used toward pixels where the 10.4 second data were 
saturated, which is only the case for some of the sources 
in NGC 6334 I. The pixel size for the images is l."2. 
Upper limits were derived from mosaics produced with 
the MOPEX program. To calculate the upper limits, 
the rms signal was found for a 5x5 pixel box centered on 

® The Submillimeter Array (SMA) is a collaborative project 
between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy &; Astrophysics of Tai- 
wan. 
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the National Science Foundation operated under agreement by the 
Associated Universities, Inc. 



the SMA sources. These values were converted to units 
of mJy per pixel, then multiplied by a factor of 5/2 to 
calculate 5cr limits for detections in a 2x2 pixel aperture. 

3. RESULTS 

3.1. NGC 6334 I 

The 1.3 mm continuum emission in NGC 6334 I is 
shown in Figure ^ We have resolved four major con- 
tinuum sources which we denote as I-SMA1..4, in de- 
scending order of peak intensity (see Table QJ. The 
continuum sources are all fairly compact. Gaussian fit- 
ting yields formal size estimates between 1.6 and 2.7 
arcsec (2800 and 4600 AU). Including compact and ex- 
tended emission, the total flux density within Figure ^ 
is 10.7 ± 0.1 Jy, equivale nt to 43 ib 5% of the single dish 
continuum flux density l|Sandel]l 1200(1) . The overall ex- 
tent of the emission (« 10") a grees well wi th the size 
estimate from single-dish maps l|Sandel]| l2000D. No other 
sources are found beyond this central concentration of 
objects. The two brightest millimeter sources (I-SMAl 
and I-SMA2) correspond within 0.5" to the two primary 
NH3 pe aks observed in the (1,1), (2,2), and (3,3| tran - 
sitions ijBeuther et alJ 120051 : iKraemer fc .T ackson"1995t). 
The H2O masers form a linear structure (Forstcr 199^ 
that appears to be associated with I-SMAl. In contrast, 
the 6.7 GHz CH3OH masers are associated with I-SMA2 
and LSMA3. LSMA3 coincides to within 0.3" of the 
peak of the 3.6 cm emissi on which tra ces a cometary 
compact Hii region (Carra ret al.11199'^ . Based on the 
3.6 cm flux density, wc estimate that the free-free contri- 
bution to the 1.3 mm flux density of I-SMA3 is ~ 1.25 Jy 
(62%). The rest of the 1.3 mm continuum emission from 
this source (and all of the other sources) originates from 
dust. LSMA3 lies 1.0" south of the NIR star IRSIE 
which has been prop osed as the exci ting source of the 
compact Hii region ijPersi et al.l 1199 6). Since free-free 
emission dominates I-SMA3, it is perhaps not surprising 
that its position lies closer to the mid-IR/cm pea k than 
to IR SIE. The know n mid-IR sources IRS- I 2 (Pcrsi et alJ 
I1998D and DPTOO 2 l|De Buizer et all 20Q2'l reside within 
the extended millimeter contours, but the lack of a com- 
pact millimet er source at the s e pos itions supports the 
hypothesis of iDe Buizer et all l)2002(l that they are not 
internally heated. With the exception of the Hii region 
nebula (I-SMA3), none of the millimeter sources are de- 
tected in the IRAC images. Indeed, no compact counter- 
part at any other wavelength has been observed toward 
LSMA4. 

3.2. NGC 6334 I(N) 

We have resolved the millimeter continuum of 
NGC 6334 I(N) into seven sources, which we denote as 
I(N)-SMA1..7 (Figure 121). Including the extended emis- 
sion, the fraction of single-dish flux recovered in our 
image is 32 ± 8%. Nearly half of the total emission 
originates from I(N)-SMA1. Although this field lacks 
any strong centimeter continuum emission, two faint 3.6 
cm (0 .3 mJy) sources have been reported (jCarral et alJ 
I2002D . One of these sources lies within 0.6" of I(N)- 
SMA4, which is also the source c l osest to a C lass II 
CH3OH maser l|Walsh et all 1199^ iNorris et al.l l!993). 
The faint continuum source with its (probable) asso- 
ciated class II CH3OH maser may represent a hyper- 
compact HII region, similar to the sources found by 
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1 7^20""54^0 53^5 53^0 

Right Ascension (J2000) 

Fig. 1. — Grayscale and black contours: 1.3mm continuum image of NGC 6334 I taken from the combined LSB plus USB data. Contour 
levels are (-3. 3. 6. 1 2. 18, 24, 36, 48, 64, 80, 100, 120, 150, 180) X 7.5 mjy/beam. Water masers are marked by the blue s ymbols: 'x' 
IMi genes et aljfl999l) and '+' ll'br stcr 1992) (shifted b y the amount sugg ested in a private communication from J. Forster in ' Carral et alj 
119971) ). Open squares mark Class II CH3OH masers (IWalsli et al.JI1998'). The black cross marks IRS IE and the six-pointed stars mark 
the three mid-infrared sources from east to west: DPTOO 2 tDe Buizer et a l. 2000), I R S-I 1 a nd IRS-I 2 Jpe Buizcr et al. 2002). The green 
crosses mark NH3 peaks (derived from images of the hyperfine lines of lBeuther et al.l 120051) '). Red contours: VLA 3.6 cm emission (-5, 5, 
15, 30, 60, 90, 120, 150 mJy beam"!). 



TABLE 1 

Properties of millimeter continuum sources in NGC 6334 I and I(N) 



Source 


J2000 coordinates 


^1 . 3mm 


-^1.3mm 




-Fs.O/jm 




a (h m s) 


5 (° ' ") 


(Jy/b) 


(Jy) 


(mJy) 


(mJy) 


I-SMAl 


17 20 53.44 


-35 46 57.9 


1.77 


3.49 ±0.70 


< 320 


a 


I-SMA2 


17 20 53.20 


-35 46 59.6 


0.96 


2.28 ±0.46 


< 240 


< 1490 


I-SMA3 


17 20 53.45 


-35 47 02.6 


0.71 


2.01 ±0.40 


a 


a 


I-SMA4 


17 20 53.12 


-35 47 03.2 


0.52 


1.10 ±0.22 


< 260 


< 1490 


Total 








10.7 ±2.1" 







I(N)-SMA1 17 20 55.21 -35 45 04.1 0.82 2.04 ± 0.41 < 0.11 < 0.71 
I(N)-SMA2 17 20 54.90 -35 45 06.8 0.35 0.50 ±0.10 < 0.35 < 0.53 
I(N)-SMA3 17 20 55.00 -35 45 07.5 0.27 0.39 ± 0.09 < 0.21 < 0.53 
I(N)-SMA4 17 20 54.69 -35 45 08.5 0.17 0.33 ± 0.07 < 0.70 < 0.77 
I(N)-SMA5 17 20 55.08 -35 45 02.0 0.27 0.28 ± 0.07 < 0.12 < 0.63 
I(N)-SMA6 17 20 54.59 -35 45 17.9 0.27 0.47 ±0.10 < 0.29 < 0.69 
I(N)-SMA7 17 20 54.96 -35 44 57.3 0.07 0.20 ± 0.05 < 0.51 < 0.68 
I(N) Total 4.6 ± 0.9" 

The extended nebula of the compact HII region precludes a meaningful infrared point source upper limit 
^ The total 1.3mm flux density includes the extended emission; the listed uncertainties include 20% calibration uncertainty 



Ivan der Tak fc Ment cn' ('2OO50 in other massive star for- 
mation regions. If this interpretation is true, then this 
object is an important power source in the I(N) region. 
Class I CH3OH masers have also been identified in this 



region at 25 GHz CBcu ther et "311120051) and at 44 GH z 
but with uncertain astrometry i Kogan fc Slvshl Il998(l . 
None of the millimeter sources have counterparts in our 
IRAC data (5(t upper limits at 4.5/im and 8/im are given 
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1 7''20'"56'0 55*0 54f0 



Right Ascension (J2000) 

Fig. 2.— 1.3 mm continuum image of NGC 6334 I(N). Contour 
levels are (-3, 3, 6, 9, 12, 15, 20, 25, 30, 35, 40, 45, 50, 60, 70, 
85, 100, 115) X 6.5 mjy/beam. Filled circles mark 44 GHz Class 
I CH3OH masers iKogan fc gTvsh 1998), the small cross marks 
a 6.7 GHz Class II CH3OH maser iNorris et al.'1993t 'Wal sh et al.l 
119981) . the X marks the faint 3.6 cm source '^arral ct al. 200 2), and 
the 6 -pointed stars mark infrared H2 knots ( Mcecath & Tieftrunk" 
[T999h . The filled triangles mark 25 GHz Class I CH3OH masers 
and the large cross indicates the NH3 peak liBeuther et al.il2005ii) . 

in Table [TJ). However, directly west of I(N)-SMA4 is 
an extended infrared source detected at 4.5 and 5.8 /im 
(FigureEJ. The elongated morphology suggests a jet or a 
reddened reflection nebulosity from an outflow cavity. In 
either scenario, the extended mid-IR emission would be 
the result of outflow activity from the embedded source 
in I(N)-SMA4. The unipolarity could be explained if 
the other lobe of the putative outflow was obscured by 
stronger extinction in its direction. This possible associ- 
ation is the only evidence, albeit indirect, for the I(N)- 
SMA sources at wavelengths < 8 /im. 

3.3. Mass estimates from the dust emission 

We can estimate the mass of the individual dust 
sources by placing our millimeter continuum measure- 
ments in the context of a si mple isotherm al model of 
optically-thin dust emission l|BeIt ran "eTaLl 12006). The 
assumption of low optical depth is justified because the 
peak continuum brightness temperature of the strongest 
dust source (TSMAl) is only 6.3 K. The input param- 
eters to the model include the dust temperature, the 
dust mass opacity coefficient at the observed wavelength 
(Ki.smm), and the gas-to-dust mass ratio. Based on our 
SMA spectra (Hunter et al. , in preparation) , we can cat- 
egorize the sources into three groups: A) strong hot-core 
line emission: I-SMAl, TSMA2 and I(N)-SMAl, B) weak 
line emission: I-SMA3 and I(N)-SMA2, and C) nearly 
line-free emission. Thus, from a qualitative viewpoint, 
the gas and dust is likely to be warmest in the group A 
sources, and coolest in the group C sources. By combin- 
ing this information with the various gas temperature 
measurements from the literature, we have assigned a 
probable value for the dust temperature of each source 
including upper and lower limits (see Table [21). 

For the sources in group A, the brightness tempera- 
ture of the arcsecond scale N H3 (2,2) emission pr ovides 
the lower-limit temperature ijBeuther et alJl2005j) . For 



TABLE 2 

Range of estimated temperatures and masses of dust 
CORES IN NGC 6334 I and I(N) 



cold limit nominal warm limit 





K 


T 


M 


T 


M 


T 


M 


Source 


(crn^ g~ 


1) (K) (Me 


) (K) (Mq) (K) 


(Mq) 


I-SMAl 


2 


75 


23 


100 


17 


300 


5.4 


I-bMA2 


2 


75 


15 


100 


11 


300 


3.5 


I-SMA3 


1 


40 


52 


60 


33 


100 


19 


I-SMA4 


1 


20 


66 


33 


36 


50 


22 


Total 






156 




97 




50 


I(N)-SMA1 


2 


65 


15 


100 


10 


300 


3.1 


I(N)-SMA2 


1 


20 


30 


40 


13 


100 


4.8 


I(N)-SMA3 


1 


20 


23 


33 


13 


50 


7.8 


I(N)-SMA4 


1 


20 


20 


33 


11 


50 


6.1 


I(N)-SMA5 


1 


20 


17 


33 


9 


50 


5.2 


I(N)-SMA6 


1 


20 


28 


33 


15 


50 


8.7 


I(N)-SMA7 


1 


20 


12 


33 


6 


50 


3.7 


I(N) Total 






145 




77 




39 



I-SMAl and I-SMA2, we set the "nominal" tempera- 
ture equal to the dust temperature (100 K) obtained 
from fits to the sin gle-dish spectr al energy distribution 
(SED) for source I fSa ndellll2000j) . We also use 100 K 
for I(N)-SMAl, because the single-dish intensity ratio 
of NH3 (6,6 ) to (3,3) indicates the presence of gas with 
T > 95 K ijKuiper et al.l 119951) : furthermore, recent in- 
terferometric measurements have localized the bulk of 
the NH3 (6,6) emission to the position of I(N)-SMAl 
(Beuther et al., in preparation). For an upper limit tem- 
perature for group A, we adopt 300 K because an exci- 
tation temperat ure of 295 K was meas ured f or H2CO in 
NGC 6334 I b vlMangum fc W ootte4 jT99l and 213 K 
for CH3OH bv lvan der TaTet al.1 Ij2003|) . For the sources 
in group B, we set the nominal tempe rature equal to 
the CO excitation temperature derived bv lKraemer et all 
l(T999|) which is 60 K for NGC 6334 I and 40 K for 
I(N). For the sources in group C, we set the nominal 
temperature equal to the value (33 K) derived from a 
Large Velocity Gra dient analysis o f a single-dish submil- 
limeter line survey ijMcCutcheon et al.ll200(^ . In both I 
and I(N), the lower limit temperature of the group B 
and C sources is taken to be the temperature of the 
coolest dust core (20 K) in the surrounding region as 
measured by jSandcU (2000), while the upper limit is 
taken to be 50 K due to the lack of direct evidence of any 
warmer gas at these positions. Measurements of Ki.smm 
in NGC 6334 I and I(N) with a 30" beam yield value s 
from 1.0 to 1.2 ± 0.6 cm g-i l|Schwartz et alJll989fl . 
However, due to the large range of temperature in the 
compact millimeter source s, we have chosen to use tabu- 
lated values of Ki.amm from lOssenkopf fc Hennin j l|1994fl 
for a density of 10^ cm~^. For sources in group A, we 
use a Ki.smm — 2 cm^ g~^ for dust grains without man- 
tles, which is appropriate for regions where protostel- 
lar heating has destroyed the solid ice but has not yet 
dispersed the dust aggregates. In the other milhmeter 
sources which show no hot core emission and are opaque 
at infrared wavelengths, the grains are likely to have thick 
ice mantles; thus we use Ki.smm = 1 cm^ g~^. Applyin g 
a gas-to-dust mass ratio (g) of 100 ijSodrosk i et al. 1997^. 
these assumptions yield the masses listed in Table|21 Due 
to uncertainties in Ki.amm and g, the uncertainties in all 
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of the masses in Table |2] are likely to be at least a factor 
of two. Subject to these uncertainties, the total mass 
estimates of the compact objects (in both the warm and 
cold temperature limits) are of similar magnitude in I 
and I(N), ranging from - 50 Mq to ~ 150 Mq. In I(N), 
this mass is divided among a larger number sources which 
are, in turn, spread over a larger angular extent. 

4. DISCUSSION 

The observed multiplicity and strength of the mil- 
limeter continuum emission is quite similar between 
NGC 6334 I and I(N), in contr ast to their strikingly 
differ ent appearance in the NIR l|Megeath fc TieftrunkI 
119991) and in ther mal hnes of NH3 and CH3OH 
(jBeuther et al.f 2005'). Although previous maser and NH3 
observations of NGC 6334 I had suggested the presence 
of multiple sources, our millimeter image provides a di- 
rect and unambiguous picture of a massive protoclus- 
ter. The estimated mass of interstellar material toward 
each source is sufficient to form a massive star and is 
significantly greater than typical objects detected in the 
single-dish surveys of low-mass star-forming regions such 
as Perseus, where an average core mass of 2.3 Mp, has 
been found with a beamsize of ^ 8000 AU ()Enoch et al.l 
12006ft . In addition, it is likely that a compact protostel- 
lar object may already be present at the center of each 
source, heating the surrounding gas and dust and lead- 
ing to our strong millimeter detections. This scenario is 
particularly likely for the strong hot-core sources I-SMAl 
and I-SMA2. It is important to note that the masses in 
Table |21 do not include the total stellar mass that may 
be present in addition to the interstellar material. 

In the case of NGC 6334 I(N), the presence of out- 
flows, masers, warm gas emission, and a faint centimeter 
continuum source also indicates the presence of embed- 
ded sources. In comparison to source I, this cluster of 
sources may be in an earlier evolutionary phase, as sug- 
gested by the lack of a NIR cluster, HII region, or bright 
mid-IR source, and the significantly high er gas column 
density inferred from N2II+ observations fPirog ov et al.l 
|2Q03). However, if the millimeter sources contain cen- 
tral protostars, then the lack of mid-IR detections to- 
ward all seven sources is especially curious, because, in 
contrast to source I, there is no confusion from bright, 
extended nebulosity in this region. If we consider the 
I(N) bolometric luminosity of 1.9 x 10"^ L0(Sandell 
120001) . the ea rliest ZAMS star t hat could be present 
is type B2V (|Hanson et alJ 11997ft . Our lowest 4.5^m 
upper limit of 0.11 mJy (from Table ^ corresponds 
to an apparent magnitude TO4.5^ni 

> 15.5. This limit 
is consistent with the non-detection of an embedded 
B2V ZAMS star as long as the extinction at 4.5/^m is 
> 4.9 mag, a ssuming an absolut e stellar magnitude of 
Mk = —0.68 l|Hanson et al.lll997ft and an intrins ic color 
{Mk - M4.5^m) = -0.08 mag fKoornneef '1983). This 
level of extinction at A.5iJ,m co rresponds to 9.3 mag of 
K band extinction (Indcbctouw et al. 2005), 83 mag of 
V band extinction (Rieke 1^ Lebofskv 19§ 5 iV and a tota l 
hydrogen column density of 1.8 x 10^^ cm~^ (|Rvteill996ft . 
Even using the smallest mass estimates from Table|3 the 
gas column density implied by the dust emission is sev- 
eral times larger than this value and can easily explain 
the non-detections. As a result, it is unclear whether 
these millimeter sources contain young protostars (B2 or 
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Fig. 3.— i.5^lm IRAC image of NGC 6334 I(N). Contours in- 
dicate the 1.3mm continuum emission (5, 10, 20, 40, 70, 100 X 
6.5 mjy/beam) and the dashed circle shows the primary beam of 
the 1.3mm SMA observations. The newly-detected infrared nebula 
happens to lie near the center of this circle, just west of I(N)-SMA4. 



later) that are simply heavily obscured, or whether they 
are intrinsically faint in the mid-IR due to their youth, 
and may eventually evolve into higher luminosity O stars. 
In either case, the high angular resolution and sensitivity 
of our observations have revealed the presence of compact 
systems where intermediate to high mass stars are being 
formed. 

In NGC 6334 I, all four of the compact millimeter 
continuum sources are concentrated within a region of 
projected radius < 0.05 pc which is signific antly smaller 
than the 0.33 pc radius of the NIR cluster ijTapia et alJ 
11996ft and the 0.23 pc radius of the search region (i.e. 
the SMA primary beam). Furthermore, the millimeter 
sources in NGC 6334 I reside close to the geometric cen- 
ter of the NIR cluster where the number density of mil- 
limeter sources (8000 pc~^) exceeds the stellar density 
of the NIR cluster (1200 pc~'^). These findings suggest 
that massive star formation is biased toward the cen- 
ter of the parent cluster. The projected separation of 
the four components of NGC 6334 I range from 5800 
to 11200 AU. Because the ratio of distances between all 
components is less than three, it satis fies the criterion 
for ide ntifying Trapezia a s described bvlAbt fc Corball^ 
()2000ft and originated bv lAmbartsumianI l)1954|) . If the 
total luminosity of NGC 6334 I (2.6 x 10^ L©) was ap- 
portioned equally to four sources, their individual lu- 
minosities ^w2uW_coiTespond to that of an 08.5 ZAMS 
star ( Han son et alJ ITQgTl) . These findings argue that 
the NGC 6334 I SMA sources may constitute a proto- 
Trapezium system, similar to the more evolved Orion 
Nebula Trapezium (10,000 AU in diameter) and many 
other optically visible nebulae and clusters. Although the 
I(N) SMA sources are more spread out, five of the sources 
of I(N) reside inside a region of 10,000 AU, suggesting 
that a similar system is forming in I(N). These obser- 
vations suggest that the presence of Trapezia in clusters 
and nebulae may result from the massive stars forming 
preferentially in the centers of clusters in the deepest part 
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of the gravitational well. 

Are t hese new observations consi s tent with previous 
claim s ijMegeath fc TieftrunkI [19991 iSollins fc MegeathI 
I2004D that NGC 6334 I(N) is in an earlier stage of 
star formation? In source I, we have identified exam- 
ples of young massive protostars in the center of an 
established infrared cluster containing low to interme- 
diate mass stars. In source I(N), we see evidence of 
massive protostars without a surrounding NIR cluster, 
which challenges the common picture of high mass stars 
forming after t he first genera tion of low mass stars (e.g. 
iHerbiel l|T96l : lKumar et alJ l|200fit) V More sensitive in- 
frared imaging is needed to probe deeper through the 
extinction for signs of lower mass protostars in this field. 
However, until such evidence is found, the question of the 
relative age of I and I(N) hinges on three main facts: 1) 
the ratio of bolometric luminosity between I and I(N) is 
large (« 140) l|Sandellll2000|l : 2) the I(N) region exhibits 
significantly less hot-core line emission than NGC 6334 I 
(Thorwirth et al. (2003); Hunter et al., in preparation); 
3) a cluster of compact dust continuum sources exists in I 
and I(N) with a comparable amount of mass (this paper). 
These observations suggest that NGC 6334 I is a more 
evolved cluster than I(N), or that I(N) is forming a clus- 
ter with a larger number of stars but of lower mass. In 
either case, our millimeter data provide strong evidence 
that NGC 6334 I(N) is forming a cluster of stars, even 
though an associated NIR cluster has not been identi- 
fied, and may be still in the process of forming. Finally, 



we note that the virial masses derived fr om single-dish 
molec ular line spectr a (including N^H"*" bvlPirogov et alJ 
ll2?)0l and HC3N bv lSollins fc MegeathI (120041) ) for both 
NGC 6334 I and I(N) are several times higher than the 
total estimated mass contained in the compact millimeter 
sources. This fact illustrates the wealth of star-forming 
material in both regions that has either not assembled 
into compact protostars or exists in a wider distribution 
of smaller cores below our sensitivity limit. It is this 
gas which may be forming the cluster surrounding the 
SMA sources in NGC 6334 I(N) and which may further 
increase the population of the cluster in NGC 6334 I. 
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